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Abstract 

A full three-dimensional angular momentum projection on top of a triaxial relativistic mean-field 
calculation has been implemented for the first time. The underlying Lagrangian is a point coupling 
model and pairing correlations are taken into account by a monopole force. This method is applied 
for the low-lying excited states in ^^Mg. Good agreement with the experimental data is found for 
the ground state properties. A minimum in the potential energy surface for the 2^ state, with 
/5 ~ 0.55,7 — 10°) is used as the basis to investigate the rotational energy spectrum as well as the 
corresponding B(E2) transition probabilities as compared to the available data. 
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One of the most discussed topics since the 60's has been the possible triaxial nature 
of 2s-ld shell nuclei. The best candidate to show such a deformation is the nucleus ^^Mg, 
which has been extensively studied both experimentally and with state-of-the-art theoretical 
models. However, the question of whether it is triaxial is still not well settled. 

Early experiments and their interpretation, such as the spectroscopic factors extracted 
from the ^'^Mg(d,p)^^Mg reaction data done by Parikh favor a prolate axial deformation. 
On the other hand, there is also evidence that supports the triaxialityof ^"^Mg, like the 
interpretation of inelastic polarized proton scattering by Lombard et al. or the existence 
of a low-lying K = 2 band However, the interpretation of this band as a rotational 
band or a band built on one-phonon 7-vibrational states, is the origin of much controversy. 
Phenomenologically, the energy spectrum for the low-lying states of ^^Mg is well described 
by the asymmetric rotor model of Davydov-Filipov with a triaxiality parameter 7 = 22°, 
while the branching ratios and the reduced transition rates are not Furthermore, the 
interband transitions are reproduced with 7 = 14°, a value that cannot give the correct 
energy spectrum [5|. 

Many microscopic studies have been done to describe the shape of ^^Mg. So far, most 
of them have relied on investigations of the potential energy surface (PES) within vari- 
ous mean-field approaches Q. Early triaxial Hartree-Fock (HF) calculations with different 
Skyrme forces showed that the shape depends strongly on the spin-orbit part of the inter- 
action [3]. With a proper treatment of this term, both Hartree-Fock-Bogoliubov (HFB) 
calculations with the Gogny force [sl as well as HF-I-BCS investigations with the Skyrme 
force jJ], yield axially symmetric prolate minima. However, a Three Dimensional (3D) 
cranked HFB calculation within the pairing-plus-quadrupole model predicts a substantial 
triaxial deformation for states with higher angular momenta J > 4 in the ground state 
band Similar investigations have been also carried out within the relativistic mean- 

field (RMF) approach, and an axial prolate ground state has been predicted for ^^Mg with 
either NLl or TMl [l^ parameter sets. 

All the mentioned microscopic studies are based on the mean-field approximation, which 
is able to take into account correlations by "spontaneous symmetry breaking", e.g., by 
the introduction of a deformed mean-field that breaks rotational symmetry in the intrinsic 
frame. In order to compare properly with the experimental data, one has to transform to 
the laboratory system and to resort to Angular Momentum Projection (AMP) to restore the 
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broken symmetry. However, due to its very demanding numerical nature, only very recently 
it has been possible to apply such a projection procedure to study the low-lying excited states 
of ^^Mg, in the context of the Generator Coordinate Method (GCM), for example with the 
Skyrme SLy4 force [3], the Gogny DIS force [l4] and the relativistic point-coupling model 
with PC-Fl force [l5|. These investigations have shown that the energy gain due to the 
restoration of rotational symmetry is of the order of several MeV. Meanwhile, it has to be 
pointed out that all these AMP+GCM studies are restricted to axially symmetric mean- 
fields. Although great success was achieved in the description of both the spectrum and 
B(E2) transition probabilities with J < 6, neither the triaxial deformation information nor 
a satisfactory description of the observed K = 2 band could be obtained. Therefore, it is 
essential to introduce triaxial quadrupole deformation at the mean-field level and perform 
full three-dimensional angular momentum projection (3DAMP). 

In the context of the phenomenological models, 3DAMP has been done using small shell 
model spaces and the corresponding effective interactions Q, H Q, Q. It is shown that 
the restoration of rotational symmetry has a strong infiuence on the topological structure 
of (/5,7) energy surface In the context of energy density functionals, 3DAMP has been 
performed on top of Hartree-Fock with a simple Skyrme-type interaction [20], or with a 
full Skyrme energy functional [2li]. In both cases cranked wave functions were projected 
to approximate a variation after projection procedure, but pairing correlations were not 
included. Only very recently, projection on particle number (PNP) and angular momentum 
with configuration mixing has been attempted in a triaxial HFB theory for the study the 
low-lying excited states of ^^Mg 22|. 

RMF theory, which relies on the ideas of effective field theory and of density functional 



theory, represents a very success: 



a few universal parameters 23, 



ul approach to low energy nuclear structure studies with 
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26!]. It incorporates many important relativistic 



effects from the beginning, such as the presence of large Lorentz scalar and vector mean 
fields with approximately equal magnitude and opposite sign that leads to a new saturation 
mechanism via the difference between the scalar and vector densities, and yields naturally 
the large spin-orbit splitting needed for the understanding of magic numbers in finite nuclei. 
In this manuscript, for the first time, 3DAMP has been implemented on top of triaxial 
relativistic mean-field theory based on point-coupling (RMF-PC) interactions, and will be 
applied in this Letter to study the low-lying excited states of ^^Mg. 
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In this approach, the intrinsic mean-field state is obtained through the minimization of 



an energy functional based on the following Lagrangian density 

—F^'F,, - ei^^7^^A^. (1) 

which contains eleven coupling constants, as, ay, cuts, oitv, Ps, Is, Iv, ^s, ^v, ^ts and 
5tv, where the subscript indicates the symmetry of the coupling: 5* stands for scalar, V for 
vector and T for isovector, while the Greek letters indicate the kind of contact interaction: 
a refers to four-fermion terms, (3 and 7 to third- and fourth-order terms respectively, and 5 
to derivative couplings. 

From the Lagrangian density in Eq.([T]), one can easily obtain an energy density func- 
tional, whose minimization provides the intrinsic mean-field state wave-function. So far, 
pairing correlations have been taken into account within the BCS approach based on a 
monopole force and a smooth cutoff factor is taken into account to simulate the effects of 
finite range [28^. To obtain the PES, the mass quadrupole moments are constrained through 
the quantities g2o and 522, which are related to the triaxial deformation parameters (3 and 7 
of the Bohr Hamiltonian by ^20 = f^-Ro/^ cos 7 and 522 = '^Rl^P sin 7, where Rq = 1.2^4^/^ 
fm. The total mass quadrupole moment q is thus given by g = A/l67r/5A/g|o + '^(122- During 
minimization, parity, D2 symmetry, and time-reversal symmetry are imposed. The densities 
are thus symmetric with respect to reflections on the x = 0, y = and z = planes. The 



parameter set chosen for the Lagrangian density in Eq.([T]) is PC-Fl [27|]. The solution of 
the equations of motion obtained from the energy density functional is accomplished by 
expansion on a set of isotropic three-dimensional harmonic oscillator basis in Cartesian co- 



ordinates with = 8 major shells. The oscillator length is chosen to be 60 = ^/h/rriLUo 



with hwo = 41 A^/^ MeV. The pairing strength constants Gn = 34.60/^4 and Gp = 33.75/^4 
are determined separately for neutrons and protons by adjusting the pairing gaps at the 
minimum of mean-field PES to the odd-even mass difference as obtained with a five-point 
formula, and kept fixed throughout the constraint calculation. 

The wavefunction I^E'q*^) with good quantum numbers J and has the form, 

K>0 

where the sum is restricted to non-negative even values of K. The index a = 1, 2, ■ ■ ■ 
labels the different collective excitations. The angular momentum projected i^-component, 
\JMK+, q), is given by 

\JMK+,q) = -^[Pi,^ + {-iyPl,_^]mq)). (3) 
where the projection operator -P/^^ is defined as {g] 

PiiK = / dnDi:iK{n)R{n), (4) 

with fl representing the set of three Euler angles. D'ljj^iVL) is the corresponding Wigner 
function. The expansion coefficients fa^{q) are determined by the generalized eigenvalue 
equation, 

E ^'^KK'iQ-, Q) - EXk'W, q)}fi''\q) = 0, (5) 

K'>0 

where the overlap kernels {O = J\f, 7i) are given by: 

Oi,K,{q]q) = {JMK + q\d\JMK' + q), d = l,H. (6) 



The overlap kernels are determined with the help of generalized Wick's theorem 29|, l30 |. 



In addition the method of Neergard and Wiist [31i is also used to determine the phase of 
the overlaps {R{Q)), and it is found that both methods are in good agreement. The norm 
overlaps for intrinsic mean-field states with deformation 7 = have been compared with 
those computed in the Gaussian overlap approximation, which turns out to be a good ap- 
proximation, as already pointed out in Ref. [15|]. This provides a very useful test of the 
numerical implementation of the angular momentum projection techniques. For the Hamil- 
tonian overlap {HR{Q)), transition (or mixed) density matrices are adopted. The numerical 
evaluation of the kernels ([6]) is carried out by a N-point Gauss-Legendre quadrature over 



the Euler angles Q = {(j),9,ip)- With Ng = 14, = = 24, it is possible to achieve a 
precision of 0.001% in the energy of a projected state with an ang ular momentum up to 
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18|, only 1/16 of the total 



J = 6. With the consideration of symmetries in the overlaps 
overlaps needs to be evaluated directly. 

The standard procedure described in Ref. ^ has been used to solve the generalized 
eigenvalue equation ([5]). By construction, the collective wave function Eq.([2]) avoids the 
appearance of zero eigenvalues, and therefore there are J/2 + 1 or (J — l)/2 collective states 
and rotational energy levels for even or odd angular momentum J 18|, respectively. These 
levels can be classified into different bands according to their B(E2) transition probabilities, 

^2 



B{E2;q,Ji,ai q,Jf,af) 



2Ji + 



J\{Jf,q\\Q2\\J^,q)\' 



(7) 



Since these quantities are calculated in full configuration space, there is no need to introduce 
effective charges, and hence e denotes the bare value of proton charge. A full detailed 
discussion about the 3DAMP+RMF-PC approach can be found in Ref. [32 1. 

Figure [1] shows the PES of ^^Mg in the deformation plane, calculated with triaxial 

RMF-PC+BCS. The minimum with E = -193.57 MeV, and = 0.48 is found which is 



less bound while compared with the data E = —198.26 MeV 
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FIG. 1: (Color online) The potential energy surface in (3-'y plane obtained by triaxial RMF- 
PC+BCS calculation for ^^Mg. The contour lines are separated by 1.0 MeV. 



The PES of the first J'" = 2"*" and J'^ = 0+ projected states for ^^Mg in the jS-'j plane 
are plotted in Fig. O There is no pronounced minimum with an obvious 7-deformation 
in the PES for the 0"*" state, which is in disagreement with the results of ref. 22|, where a 
calculation with 3DAMP+PNP based on non-relativistic Skyrme HFB shows a pronounced 
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FIG. 2: (Color online) Three-dimensional angular momentum projected energy surfaces with J = 
(left panel) and J = 2 (right panel) in the /3-7 plane for ^''Mg. The contour lines are separated by 
1.0 MeV 
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FIG. 3: Comparison between the theoretical and experimental bands in ^^Mg. The excitation 
energies are normalized to the ground state and given in MeV. The widths of the arrows represent 



the B{E2 I) values. The experimental data is taken from Refs. js] and 
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triaxial minimum with (3 = 0.6 and 7 ^ 16° for the O"*" state. The energy gain due to the 
restoration of rotational symmetry in the ground state of ^^Mg in the present investigation is 
about 5 MeV. The minimum on the PES for the 0"^ state has a larger deformation jSp = 0.56 
close to the data P = 0.57 |3J]. Keeping in mind the strong pairing gaps found in our mean- 
field calculations, the introduction of PNP in our calculations is not expected to introduce 
great changes in the results. 

The right panel of Fig. [2] shows that there is a minimum in the PES for the 2^^ state 
localized at /5 ^ 0.55, 7 ^ 10°. Using the 3DAMP+RMF-PC approach we calculate the 
spectrum of ^^Mg based on this triaxial mean-field state. This is the only way to obtain 
K = 2 bands in our calculation. In Fig. [3] we show this spectrum and compare it with 
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the known experimental data. The observed energies are quahtatively reproduced by the 
theory, even though the calculated levels are systematically higher in energy, which might 
be due to the fact that all the excited states are projected from the same mean-field state. 
The projected energy spectrum would be much more compressed with the inclusion of a 
cranking term in the mean-field calculation 
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20, |2l|]. This is a strong indication that, for 



a quantitative understanding of the spectra, the variation should be carried out after the 
projection. Of course this goes beyond the present computer capabilities. 

In conclusion, we have studied the properties of the ground state and low-lying excited 
states in ^^Mg using a full three dimensional angular momentum projection on top of a 
relativistic mean-field calculation based on a point-coupling Lagrangian (3DAMP-I-RMF- 
PC) with monopole pairing. Ground state properties are found to be reproduced rather well 
with this approach, even though no pronounced minimum with obvious triaxial deformation 
is present. A minimum with (3 ~ 0.55, 7 ~ 10° has been found on the PES of the first 
2^ state, on top of which the experimentally observed excitation energies and the B(E2) 
transition probabilities can be qualitatively reproduced. However, the predicted spacing 
between the levels is overestimated by this approach. It is expected that this can be cured 
by the introduction of a cranking term at the mean-field level. These results hint that a 
full GCM approach is needed for a full understanding of the properties of low-lying states 
in ^^Mg, and work towards this goal is in progress. 
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